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Abstract 

The potential energy surfaces of even-even 146_156 Sm are investigated in the constrained 
reflection- asymmetric relativistic mean- field approach with parameter set PK1. It is shown that 
the critical-point candidate nucleus 152 Sm marks the shape/phase transition not only from U(5) to 
SU(3) symmetry, but also from the octupole-deformed ground state in 150 Sm to the quadrupole- 
deformed ground state in 154 Sm. By including the octupole degree of freedom, an energy gap 
near the Fermi surface for single-particle levels in 152 Sm with f$2 = 0.14 ~ 0.26 is found, and the 
important role of the octupole deformation driving pair f2f 7 / 2 and ^ lz 13 / 2 is demonstrated. 

PACS numbers: 21.10.-k, 21.60.Jz, 27.70. +q 



* Email: sqzhang@pku.edu.cn 
< Email: mengj@pku.edu.cn 



1 



I. INTRODUCTION 



The first-order phase transition between spherical U(5) and axially deformed SU(3) 
shapes has received widespread attention in the past decade. It was shown that Sm 
and other N = 90 isotones are the empirical examples of the analytic description of nuclei 
at the critical point of such a transition yl . Theoretical studies on the phase transition were 

n 

typically based on phenomenological geometric models of nuclear shapes and potentials [2J 

and algebraic models of nuclear structure |4j]. The first calculations, establishing a link 

between dynamical symmetry models and microscopic theories, were carried out using the 

relativistic mean-field (RMF) approximation in the Sm isotopes [5]. Along this line, much 

work was done using either relativistic or nonrelativistic approaches 10-12]. 

Normally the regions of nuclei with strong octupole correlations correspond to either 

the proton or neutron numbers close to 34 (lgg/2 2p 3 / 2 coupling), 56 (I/Lli/2 2c?5/2 

I I 

coupling), 88 (1*13/2 -H- 2/ 7 / 2 coupling), and 134 (lj'15/2 -H- 2g 9 / 2 coupling) [13]. A variety 
of approaches were applied to investigate the role of the octupole degree of freedom in Sm 
and the neighboring nuclear region. The Woods-Saxon-Bogoliubov cranking model is used 
to study the shapes of rotating Xe, Ba, Ce, Nd, and Sm nuclei with A^ = 84 — 94 and the 
expectations of octupole-deformed mean fields at low and medium spins are confirmed [ijj . 
The spdf interacting boson model is applied to Sm isotopes with A^ = 86 — 92 to examine the 
signatures of octupole correlations [15]. Based on a collective rotation- vibration Hamiltonian 
in which the axial quadrupole and octupole degrees of freedom are coupled, the energy levels 
and electromagnetic transition probabilities for A^ = 90 isotones are well reproduced [l^ . 

Very recently, a new K n = 0~ octupole excitation band was observed in 152 Sm and a 
pattern of repeating excitations built on the 0^ level similar to those built on the ground 
state emerges \x^. It was suggested that 152 Sm, rather than a critical-point nucleus, is a 
complex example of shape coexistence [l7| • 

Based on the investigations mentioned previously, it is timely and necessary to investigate 
the Sm isotopes in a microscopic and self-consistent approach with the octupole degree of 
freedom. The newly developed reflection-asymmetric relativistic mean-field (RAS-RMF) 
approach is a good candidate for this purpose [3] considering the remarkable success of 



RMF theory 
exotic nuclei 



22, 



in describing many nuclear phenomena related to stable nuclei 



23], as well as supernova and neutron stars 24]. In Ref. |l8j], the RAS- 
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RMF approach was first applied to the well-known octupole-deformed nucleus 226 Ra, and 
reproduced well both the binding energy and deformation. 

In this article, the RAS-RMF approach will be applied to investigate the potential energy 
surfaces (PES) of even-even 146 - 156 g m isotopes in the (@2, $3) plane and the shape evolution 
involving the octupole degrees of freedom will be analyzed. 



II. FORMALISM 



The basic ansatz of the RMF theory is a Lagrangian density where nucleons are described 
as Dirac particles that interact via the exchange of various mesons and the photon. The 
mesons considered are the isoscalar-scalar cr, the isoscalar- vector u, and the isovector- vector 
p. The effective Lagrangian density reads [25] 



■0 



iYdf, - M - g a o - g^u^ - g p ^T ■ p M - ey 
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(1) 



in which the field tensors for the vector mesons and the photon are, respectively, defined as 







- duUp, 




= d^Pv 


- dvp^ 


F 

- 1 flV 


= d,A u 


- d u A, 



(2) 



Using the classical variational principle, one can obtain the Dirac equation for the nucle- 
ons and the Klein-Gordon equations for the mesons. To solve these equations, we employ 
the basis expansion method, which was widely used in both the nonrelativistic and relativis- 
tic mean-field models. For axial-symmetric reflection-asymmetric systems, where nonaxial 
quadrupole and octupole deformations are plainly excluded, the spinors are expanded in 
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terms of the eigenfunctions of the two-center harmonic-oscillator (TCHO) potential 

/ 1 



V(r ± ,z) = ^Muiri+l 



-Mu{(z + zi) 2 , z < 



-Mu 2 (z - z 2 )\ z^O 



(3) 



where M is the nucleon mass, z\ and z 2 (real, positive) represent the distances between 
the centers of the spheroids and their intersection planes, and are the corresponding 

oscillator frequencies for z < (z ^ 0) 18]. The TCHO basis was widely used in stud- 
ies of fission, fusion, heavy- ion emission, and various cluster phenomena 26]. By setting 
proper asymmetric parameters, the major and the z-axis quantum numbers are real num- 
bers very close to integers, and the integers are used in the Nilsson-like notation Q[Nn z mi] 



18|- 



for convenience. More details can be found in Ref. 

The binding energy at a certain deformation is obtained by constraining the mass 
quadrupole moment (Q 2 ) to a given value \i% (27| 



(H) + -C«Q 2 > - 



(4) 



where C is the curvature constant parameter and fi 2 is the given quadrupole moment. 

The expectation value of Q 2 is (Q2) = (Q2)n + {Qdp witn {Qz)n& = (2r 2 P 2 (cos 9)) n}P . 

3 

The deformation parameter (3 2 is related to (Q 2 ) by (Q 2 ) = —^=Ar 2 (3 2 with r = i^o^l 1 ^ 3 



V57T 

(Ro = 1.2 fm) and A the mass number. The octupole moment constraint can also be applied 

3 

similarly with (Q 3 ) = (Q 3 )n + (Q3) P , {Qz)n, P = (2r 3 P 3 (cos#)) n!P , and (Q 3 ) = ^=Ar 3 p 3 . By 
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constraining the quadruple moment and octupole moment simultaneously, the total energy 
surface in the (/3 2 , ^3) plane can be obtained. 



III. RESULTS AND DISCUSSION 



The properties of even-even 
proach with parameter set PK1 



146-156 



Sm are calculated in the constrained RAS-RMF ap- 



The parameter set PK1 is one of the best parameter 
sets available in the framework of RMF theory, which, as usual, is obtained by fitting the 
masses of selected spherical nuclei as well as saturation properties of nuclear matter. The 
success of universal RMF parameter sets was demonstrated for describing the properties of 
spherical [l9| and deformed nuclei [20|, |2l| , and they are believed to be appropriate for appli- 



cation in octupole-deformed nuclei. In return, the application for octupole-deformed nuclei 



TABLE I: The total binding energy (in MeV) as well as the quadrupole deformation 02 and 
octupole deformation {3$ of the ground states of even-even 146_156 Sm obtained in the constrained 
RAS-RMF approach with PK1, in comparison with the available experimental data. 



Nucleus 


^cal 


ocal 
P2 


ocal 
P3 


E ex P [29J 


/3 2 exP [30] 


146 Sm 


1213.38 


0.09 


0.07 


1210.91 




148 Sm 


1227.27 


0.14 


0.08 


1225.39 


0.14 


150 Sm 


1241.59 


0.18 


0.13 


1239.25 


0.19 


152 Sm 


1254.46 


0.20 


0.15 


1253.10 


0.31 


154 Sm 


1267.76 


0.32 


0.00 


1266.94 


0.34 


156 Sm 


1280.08 


0.33 


0.02 


1279.99 





will also provide a further test. The TCHO basis with 16 major shells for both fermions 
and bosons is used. The pairing correlations are treated by the BCS approximation with a 
constant pairing gap A = 11.2/y/A MeV. 

The binding energy, quadrupole, and octupole deformations are listed for the ground 
states of 146_156 Sm in Table HI The binding energies are well reproduced within 0.2%. More- 
over, excellent agreement is obtained for the quadrupole deformations except 152 Sm (see the 
discussion in the following). 

To investigate the shape evolution in the Sm isotopes, the contour plots of the total 
energies as functions of 02 and (3% for 146_156 Sm are shown in Fig. [H which have up-down 
symmetry in the (/?2, fiz) plane because of the equivalence between the states with positive 
and negative fi^. 

It is shown in Fig. [T] that for the ground states, 146 ' 148 Sm are near spherical, 150 Sm oc- 
tupole deformed, and 154,156 Sm well deformed while 152 Sm marks the transition from octupole 
to quadrupole deformed. In detail, for 146 Sm, with iV = 84 close to the magic number 
82, the ground state is near spherical with (/3 2 , / 5 3 ) = (0.08, 0.08). For 148 Sm and 150 Sm, 
with increasing neutron number, the deformations (fa, P3) gradually increase. Particu- 
larly for 150 Sm with iV = 88, a global minimum with substantial quadrupole and octupole 
deformations^, /?3) = (0.19, 0.14) is well developed, which is about 1.36 MeV deeper than 
the corresponding quadrupole-deformed state. For 152 Sm, the global minimum moves to 
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FIG. 1: (Color online) The contour plots of total energies for the even-even 146 ~ 156 Sm in (/?2, fy) 
plane obtained in the RAS-RMF approach with PK1 and constant- A pairing. The energy separa- 
tion between contour lines is 0.25 MeV. The global minima and other local minima are denoted by 
"•" and respectively. 



(P2, /#3) = (0.20, 0.15), while a quadrupole minimum emerges at (/?2, /?3) = (0.29, 0). It is 
interesting to note that the deformation of this minimum is quite close to the experimental 



value /?2 Xp 



0.31 



30j listed in Table U The energy difference between the two minima 



is 0.33 MeV with a 0.5 MeV barrier in between. For 154 ' 156 Sm, the ground states are well 
quadrupole deformed with (/3 2 , ^3) ~ (0.33, 0). Similar PES's can also be obtained with 
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other parameter sets such as NL3 31[. In addition, one notes that the oblate minima shown 
in Fig. [T]may not be stable against the 7 direction {9]. 

In Ref. 5j, axially deformed RMF calculation with a variety of effective interactions was 
performed for 144_158 Sm to discuss the transition from spherical U(5) to axially deformed 
SU(3) shapes. It was shown that the PES's of 144 ' 146 Sm are minimized near spherical and 
of 154_158 Sm well-deformed, while in between the PES's of 148 > 150 > 152 Sm are found to be 
relatively flat. From Fig. [H it is shown that the conclusion remains to be true even with the 
inclusion of the octupole degree of freedom. 

Meanwhile, we do find something new for 152 Sm after including the octupole degree of 
freedom, that is, 152 Sm marks the shape/phase transition not only from U(5) to SU(3) 
symmetry, but also from the octupole-deformed to quadrupole-deformed case. 

Quite recently in Ref. 17J, a new = 0~ band was observed, which has a remarkable 



similarity in its El transition to the first excited K 71 = + band as the lowest K 71 = 0~ band 
to the ground-state band. A pattern of repeating excitations built on the level similar 
to those built on the ground state is claimed to indicate that 152 Sm is a complex example 
of shape coexistence rather than a critical-point nucleus. These observations can be well 
understood from the PES obtained previously. For 152 Sm with the octupole minimum at 
/33)=(0.20, 0.15) and the quadrupole minimum at (/?2, /33)=(0.29, 0), if one performs 



the generator coordinate method (GCM) 27J, |32J calculation with the PES, two low-lying 



states in the {^fis) plane with similar quadrupole deformation will be obtained, which are a 
mixture of quadrupole and octupole deformation configurations. Based on these two states, 
the pattern of repeating excitations is expected. 

To understand the evolution of the octupole deformation microscopically, the neutron 
single-particle levels in 152 Sm for the states minimized with respect to /?3 and the states 
with /3 3 = for /3 2 = 0.14 ~ 0.26 are shown in Fig. [2j The levels near the Fermi surface 
are labeled by Nilsson-like notations Q[Nn z mi] of the largest component. In the left panel 
of Fig. El a large energy gap with N = 88 near the Fermi surface is found, which is related 
to the softness of the potential energy surface in the quadrupole and octupole degrees of 
freedom in 152 Sm. There is no obvious neutron gap near the Fermi surface for the states 
with /3 3 = 0. In addition, the proton single-particle levels for the corresponding states are 
shown in Fig. [31 where no obvious energy gaps can be found near the Fermi surfaces. 

It is well known that for nuclei with N ~ 88 or Z ~ 56 the octupole deformation 



7 



0.14 0.16 0.18 0.20 0.22 0.24 0.26 
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P 2 



FIG. 2: (Color online) Neutron single-particle levels of 152 Sm in RAS-RMF approach with PK1 as 
functions of /?2 for states minimized with respect to /3% (left panel) and states with = (right 
panel). The dash-dot lines denote the corresponding Fermi surfaces. The levels near the Fermi 
surface are labeled by Nilsson-like notations n[Nn z mi] of the first component at fii = 0.20. The 
corresponding /?3 are shown in the inset. 

driving pairs of orbitals include (1/2/7/2, ^1*13/2) and (Tr2d 5 / 2 , 7rl/i n / 2 ), which in the ax- 
ially deformed case will be subgrouped as (z/l/2[541], i/l/2[660]), (i/3/2[532], i/3/2[651]), 
(i/5/2[523],i/5/2[642]), (v7/2[5U), z/7/2[633]), and (tt1/2[431], tt1/2[550]), (tt3/2[422], 
7r3/2[541]), (7r5/2[413], 7r5/2[532]), respectively. It is interesting to investigate the per- 
formance of such pairs in the single-particle levels near the Fermi surfaces in Figs. [2] and [3j 
These levels together with their BCS occupation probabilities and corresponding contribu- 
tions from the four leading components are shown in Table HT1 Taking the level z/3/2[521] as 
an example, its second (20.1%) and third (15.5%) components compose an octupole deforma- 
tion driving pair (i/3/2[532], i/3/2[651]). Similarly, one can find the pair (5/2 [523], 5/2 [642]) 
for i/5/2[523], the pair (3/2[532], 3/2[651]) for i/3/2[532], and the pair (1/2[541], 1/2 [660]) 
(the fifth component 1/2 [660] with 6.0%, which is not listed in Table [XT]) for i/l/2[530]. 
However, for the proton side, no octupole deformation driving pairs are found among the 
four leading components. Therefore the neutron orbital takes a more important role in the 
evolution of octupole deformation in Sm isotopes than the proton one, consistent with the 



S 



energy gaps presented in Fig. [2j 




IV. CONCLUSION 

In this article, the PES's of even-even 146_156 Sm in the (fa, ^3) plane are obtained by the 
constrained RAS-RMF approach, and the single-particle levels near the Fermi surfaces for 
the nucleus 152 Sm are studied. It is shown that the critical-point candidate nucleus 152 Sm 
marks the shape/phase transition not only from U(5) to SU(3) symmetry, but also from the 
octupole-deformed ground state in 150 Sm to quadrupole-deformed ground state in 154 Sm. 

Furthermore, the microscopic PES for the nucleus 152 Sm is consistent with the claimed 
shape coexistence based on the observation of repeating excitations built on the + level 
similar to those built on the ground state [171 ]. 

By including the octupole degree of freedom, an energy gap near the Fermi surface for 
single-particle levels in 152 Sm with N = 88 and fa — 0.14 ~ 0.26 is found, which is related 
to the softness of its potential energy surface in the quadrupole and octupole degrees of 
freedom. From the energy gap and the components of the single-particle levels near the 
Fermi surface, it is demonstrated that the neutrons play an important role for the octupole 
deformation driving in 152 Sm. 
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TABLE II: Single-particle levels near the Fermi surface for the ground state with (/?2, fa)= (0.20, 
015) in 152 Sm together with their BCS occupation probabilities and corresponding contributions 
from the four leading components. The components originating from the octupole deformation 
driving pairs of orbitals (v2f 7 / 2 , flii.3/2) an d (^2^5/2, ^1^11/2) are m bold. 



level 


occu. 


1st comp. 


2nd comp. 


3rd comp. 


4th comp. 


i/3/2[521] 


0.401 


3/2[521] 


33.3% 


3/2 [532] 


20.1% 


3/2[651] 


15.5% 


3/2[631] 


9.4% 


z/5/2[523] 


0.434 


5/2 [523] 


57.2% 


5/2 [532] 


16.0% 


5/2 [642] 


6.7% 


5/2[633] 


5.2% 


z/3/2[532] 


0.946 


3/2 [532] 


46.4% 


3/2[541] 


21.1% 


3/2[512] 


8.8% 


3/2[651] 


5.8% 


i/l/2 [530] 


0.947 


1/2 [530] 


33.9% 


1/2[541] 


18.1% 


1/2[510] 


8.0% 


1/2[651] 


6.2% 


vrl/2[541] 


0.219 


1/2 [420] 


21.4% 


1/2[541] 


21.4% 


l/2[440] 


17.7% 


1/2[521] 


9.8% 


vr7/2[404] 


0.763 


7/2 [404] 


85.9% 


7/2[413] 


7.3% 


7/2 [514] 


2.8% 


7/2 [604] 


0.4% 


vr3/2[541] 


0.834 


3/2[541] 


34.6% 


3/2[411] 


20.0% 


3/2[521] 


17.9% 


3/2[532] 


9.1% 


vr5/2[413] 


0.952 


5/2[413] 


76.0% 


5/2 [422] 


11.2% 


5/2[523] 


4.7% 


5/2 [303] 


2.0% 
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